A New Nucleosynthesis Constraint on the Variation of G 
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Big Bang Nucleosynthesis can provide, via constraints on the expansion rate at that time, limits 
on possible variations in Newton's Constant, G. The original analyses were performed before an 
independent measurement of the baryon-to-photon ratio from the cosmic microwave background 
was available. Combining this with recent measurements of the primordial deuterium abundance in 
quasar absorption systems now allows one to derive a new tighter constraint on G without recourse 
to considerations of helium or lithium abundances. We find that, compared to todays value, Go, 
Gbbn /Go = Lollo p at the 68.3% confidence level. If we assume a monotonic power law time 
dependence, G oc t~ a , then the constraint on the index is —0.004 < a < 0.005. This would translate 



into 



-3 x 10" 18 yr- < (G/G) today 



< 4 x 10" 



The predictions of the light element abundances from 
big bang nucleosynthesis (BBN) has long served as a 
powerful probe of the early Universe. These predictions 
depend on a number of parameters such as the num- 
ber of light particle species in equilibrium at that time, 
the baryon to photon ratio (77), nuclear reaction rates, 
and the gravitational constant, G. BBN constraints on 
cosmological parameters have been traditionally limited 
by uncertainties in the observed primordial abundances 
of the light elements as well as degeneracy among the 
independent parameters. In particular, without a tight 
independent constraint on the baryon density, variation 
in the predicted abundance of light elements could be 
largely compensated, at least within observational un- 
certainties, by varying 77. 

Recently, WMAP has provided an independent mea- 
sure of 77 U allowing one to use the baryon density as 
an independent constraint on BBN instead of inferring 
its value using the apparent agreement between light el- 
ement abundance predictions and observations. In this 
vein, prior to the WMAP results other CMB experiments 
had sufficiently constrained 77 to allow tests for the con- 
sistency of BBN |3| . 

At the same time recent measurements of the primor- 
dial deuterium abundance (D/H) using quasar absorp- 
tion by intervening high redshift systems has dramati- 
cally altered its significance in comparing BBN predic- 
tions with observations Q. 

Combining these two new observations gives one a 
powerful new handle to use in order to constrain cos- 
mological parameters using BBN predictions. Here we 
explore the impact of these developments on one's ability 
to constrain the variability of G. Variation of the grav- 
itational constant was originally postulated by Dirac Q 
and remains a key component of many theories that seek 
to resolve various hierarchy problems in particle physics. 
Varying the gravitational constant has significant impact 
during BBN. 

While previous analyses (for example jHQ) attempted 



to utilize all BBN abundance predictions as a way to limit 
the effect of uncertainties in 77, as can be seen from the 
results of , there is good reason to believe that utiliza- 
tion of deuterium alone, now that its primordial value 
has been tightly constrained, might be sufficient to yield 
a stronger constraint. This is because deuterium produc- 
tion is particularly sensitive to the value of the expansion 
rate during BBN, which in turn depends upon the value 
of G. While helium is also sensitive to this rate, existing 
systematic uncertainties in its primordial value confuse 
the situation. In addition, while earlier analyses were 
performed before it was empirically demonstrated that 3 
light neutrino species exist, this residual uncertainty has 
now been removed. 

It has been known for almost three decades that deu- 
terium is only produced in significant quantities in the big 
bang This, coupled with the fact that deuterium 

depends sensitively on 77, makes it an excellent probe of 
the Universe at the time of BBN (for example see 0). 
The detection of deuterium in a number of high redshift 
quasar absorption systems allows for a precise determi- 
nation of the primordial deuterium abundance. Kirkman 
et al. 1QJ give as the best current estimate of the pri- 
mordial deuterium abundance D/H — 2.78lg'gg x 1CP 5 
based on the observation of five quasar absorption sys- 
tems. This corresponds to log (D/H) = -4.556 ± 0.064 
where the errors are assumed to be Gaussian. 

The new CMB constraints on 77 are derived from 
the temperature of the CMB and a fit of the power 
spectrum to a set of parameters defining a cosmologi- 
cal model. Here we are interested in fl^h 2 where h = 
-f/o/100 kms" 1 Mpc -1 and ilb = pb/Pc> the ratio of the 
baryon density to the critical density. A combined anal- 
ysis including WMAP, CBI, ACBAR, 2dFGRS, and Ly- 
man a data gives Q b h 2 = 0.0224 ± 0.0009 for ACDM 
models for the two choices of either a power law spec- 
tral index or a running spectral index 0. Combined 
with the photon density, n 7 = 410.4 ± 0.9 cm -3 , as de- 
termined from the COBE temperature measurement 
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FIG. 1: Observational limits and theoretical expections for 
D/H versus rj. The one (light shading) and 2 (dark shad- 
ing) sigma observational uncertainties for D/H and r\ are 
shown. They do not appear as ellipses due to the linear scale 
in D/H but logarithmic uncertainties from the observations. 
The BBN predictions are shown as the solid curves where 
the width is the 3% theoretical uncertainties. Three different 
values of Gbbn /Go are shown. 

gives rj — (6.13 ±0.25) x 1CP 10 where we again assume 
Gaussian errors. 

To determine the deuterium abundance as a function 
of rj and G from BBN we use the standard Kawano code 
[l2|. We use updated reaction rates [ill as well as used 
the latest value for the neutron half-life |b| . Krauss and 
Romanelli first demonstrated the need to incorpo- 
rate reaction rate uncertainties using Monte Carlo tech- 
niques if one is to properly derive BBN constraints in 
general. The residual uncertainty in the predicted deu- 
terium abundance is quite small. We have included this 
theoretical uncertainty in in the likelihood analysis as an 
independent Gaussian with a constant 3% uncertainty at 
the one sigma level [IE E3l That is > ^bbn = 0.03(£>/#) 
where (D/H)(t],G/Gq) comes from the Kawano code 
as noted. The uncertainties in G today are sufficiently 
small 0| and have not been included. In addition, due 
to the short time interval associated with BBN we have 
assumed that G remained constant throughout the pe- 
riod of light element production. We have verified that 
our results are unchanged if we allow G to vary during 
this time according to G oc t~ a . Figure ^ shows the 1 
and 2 sigma ellipses for the observed values of D /H and 
rj. Also shown are curves for the BBN predictions for 
(D/H)(r),G/G ) for three values of G/G with the 3% 
theoretical uncertainties from Monte Carlos of the reac- 
tion rates. 

With the assumption that the observational and the- 
oretical errors follow Gaussian distributions as dis- 
cussed above, we assign a joint likelihood function 
C(D/H,t],G/Gq) where D/H is a function of G and 
rj. The likelihood distribution for G (which can be nor- 



malized to give the probability distribution) is found by 
marginalizing over D/H and rj. Applying this procedure 
we find G/G a = l-01±g;fg at the 68.3% confidence level, 
G/G = l-Oltolo at the 95% confidence level. 

Assuming a monotonic power law time dependence 
G oc t~ a one finds the power law index —0.004 < a < 
0.005 at the 68.3% confidence level, -0.009 < a < 0.010 
at the 95% confidence level. In this case one infers 

-3 x 10~ 13 yr- 1 ( G/g) < 4 x 10~ 13 yr" 1 . This 

V / today 

is over an order of magnitude stronger than direct con- 
straints that can be obtained on the variation of G today. 

Our new constraint is about a factor of two stronger 
than the previous BBN constraint on G/G [||, which is 
significant, given that we utilize only deuterium in our 
analysis, and do not exploit other light element abun- 
dance predictions. In particular, as observational con- 
straints on these elements continue to improve, along 
with further improvements on the deuterium abundance 
measurements and the baryon to photon limit one can 
expect BBN constraints on variability of G to continue 
to improve significantly. 

LMK thanks G. Steigman for useful discussions, and 
pointing out early work on this subject. 
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